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The Source Galerkin method finds approximate solutions to the functional differential equations of field theories 
in the presence of external sources. While developing this process, it was recognized that approximations of the 
spectral representations of the Green's functions by Sine function expansions are an extremely powerful calculative 
tool. Specifically, this understanding makes it not only possible to apply the Source Galerkin method to higher 
dimensional field theories, but also leads to a new approach to perturbation theory calculations in scalar and 
fermionic field theories. This report summarizes the methodologies for solving quantum field theories with the 
Source Galerkin method and for performing perturbation theory calculations using Sine approximations. 



1. SOURCE GALERKIN METHOD 

To analyze a theory, we start with a Lagrangian in 
the presence of external sources. The functional 
differential equations satisfied by Z are generated 
in the usual way. The Source Galerkin technique 
is a systematic, iterative algorithm to solve these 
equations to increasing accuracy. The lowest or- 
der solution is given by an ansatz which is an 
expansion over the external sources of the theory, 
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We determine parameters of our expansion by im- 
posing the constraints, 
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where the number of constraints matches the 
number of free parameters. If Z* is equal to the 
exact generating functional, Z, then these condi- 
tions are always true. 

The integrals we encounter are of the form 
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We choose G xy and H^xyz so that the resid- 
ual associated with the equation(s) of motion, 
EZ* — R(j), is as small as possible. We 
do this systematically by using weighted resid- 
uals 0. We introduce a set of test functions 
{j(x),j(x)j(y),j(x)j(y)j(z), . . . } and the mea- 
sure 
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Fermions are dealt with in essentially the same 
way, using anticommuting sources This 
method can be used in two or fewer dimensions, 
but memory requirements become outrageous in 
higher dimensions J3|||. To surmount this obsta- 
cle, the problem needs to be broken into smaller 
pieces and we must utilize the symmetries of the 
action. We do this by rewriting the usual two 
point Green's function, observing that any exact 
two point function can be represented as a sum 
of free two point functions: 
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Introducing A as a cutoff, these are represented 
as 

1 



Gbcc{x) = J^yj dp -p 



Using Sine functions ||||, we can write this in a 
more numerically tractable form, 
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Reducing the free propagator to sums over Gaus- 
sians suggests a parameterization of Z* . We write 
In Z* as a Feynman sum of graphs with external 
legs replaced by sources but, unlike a perturba- 
tion theory expansion, the free propagators have 
arbitrary mass and normalization. The method 
has been applied to various field theories on the 
lattice and on the continuum. The ansatz is 
Poincare invariant and will obey other symme- 
tries of the theory by construction, giving ample 
freedom to span the solution space. The proce- 
dure can be summarized as follows: 

1. Start with an ansatz solution Z* and derive 
the Source Galerkin equations by inserting 
the ansatz into the Schwinger-Dyson equa- 
tions, and set the parameters. This involves 
functional and space-time integrals and pro- 
duces nonlinear equations. 

2. To improve the ansatz, expand Z* with 
more terms and "iterate" . 

A full account of the Source Galerkin approach 
is currently in preparation. We have made sig- 
nificant progress beyond the initial lattice imple- 
mentation of Source Galerkin method @J§]. We 
now can calculate on the continuum. Applica- 
tions to the <f> 4 in D = 1 and 2, and four fermion 
calculations |],|| give excellent results, even with 
simple Z*'s which include only two terms. The 
two-dimensional nonlinear cr-model Q results are 
compared to 1 /N calculations and the leading or- 
der (5 function calculations are identical. More- 
over, treating this model provided a laboratory 



for investigating renormalization issues in Source 
Galerkin calculations. A study of the symmetry- 
breaking sector of Higgs model is planned for the 
future. 

2. PERTURBATION THEORY 

The Sine function representation used above also 
leads to a new method for numerical evaluation 
of Feynman diagrams in conventional perturba- 
tion theory. In turn, the perturbative analysis is 
a powerful guide to "exact" Source Galerkin cal- 
culations. This approach uses generalized Sine 
functions to approximate the propagators as in- 
finite sums whose spatial or momentum depen- 
dence appears only in terms like exp(— x 2 ). Con- 
sequently, all integrations over internal momenta 
or vertex locations are reduced to Gaussian in- 
tegrals, which can be performed analytically. To 
compute the diagram one has to numerically eval- 
uate a rapidly converging sum with dimension- 
ality equal to the number of propagators in the 
graph. The Sine function representation has been 
applied to scalar and fermionic field theories in 
P~pT[| . Here we review the application of this ap- 
proach to QED calculations. 

As shown in fll|| , the fermion propagator is ap- 
proximated as 
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while the Feynman gauge photon propagator is 
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The accuracy and speed of any computation per- 
formed using these expressions are governed by 
three factors. First, some error is introduced 
through the Sine function method itself. This er- 
ror decreases exponentially with h and for h — * 
the approximation becomes exact. Additional 
error comes from the truncation of the infinite 
sums in the approximate versions of the propaga- 
tors. The general term of the sum, and thus the 
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Table 1 

Correction to the magnetic moment calculated 
with different values of h. 



h 


Correction xck/tt 


Error xa/ir 


0.2 


0.4999999999949221 


5.1 x 10"^ 


0.4 


0.5000000000008583 


-8.6 x 10~ 13 


0.6 


0.500000008001 


-8.0 x 10~ y 


0.8 


0.499999900285 


9.97 x 10-* 


1.0 


0.499998 


2.3 x 10"° 


1.2 


0.49977 


0.00023 


1.4 


0.4987 


0.0012 



truncation error, decreases exponentialy with the 
summation parameter k. Finally, finite numerical 
precision limits the maximum accuracy achiev- 
able in any computation. Setting a value of h fixes 
the expected precision of the result. Since both 
errors due to truncation of sums and the Sine ex- 
pansion decrease exponentially with correspond- 
ing parameters, the accuracy of the whole com- 
putation improves exponentially as h decreases 
while amount of work necessary to perform the 
summations with required precision increases lin- 
early. 

In order to test this method we have evaluated 
several second and fourth order diagrams. The 
results of computation of second order contribu- 
tion to the anomalous magnetic moment of an 
electron are shown in Table |[ This illustrates 
that accuracy as high as one part in 10 13 can be 
achieved using the method in question on 64-bit 
hardware, and this limit dominates for h < 0.3. 
Performing the same calculation with 32-bit ac- 
curacy yields a maximum accuracy of one part in 
10 6 , for h = 0.6. 

Our initial tests show that the Sine representa- 
tion method is capable of producing results whose 
accuracy is limited only by numerical precision of 
the hardware in a fraction of the time required 
to obtain matching results by more conventional 
methods. 



3. SUMMARY 

The Source Galerkin method promises to tackle 
a range of problems which are difficult or impossi- 
ble to formulate using conventional Monte Carlo 
approaches. Our methods work on the continuum 
and avoid most of the usual fermionic problems. 
They are well suited for calculation in any phase 
of a theory and are not susceptible to the "sign 
problem" . 

Sine function representations can be used as 
a universal tool for performing conventional per- 
turbation theory calculations. Accuracy as high 
as one part in 10 15 was achieved during testing. 
The method is extremely fast and can be used to 
quickly obtain estimates of the final result. High 
precision computations done using the Sine repre- 
sentation take significantly less time than compa- 
rable calculations performed by with Monte Carlo 
methods. We are working on the tasks of au- 
tomating this method and applying it to higher 
order diagrams, as well as higher iterations of 
Source Galerkin. 

Computational work was performed at the 
Theoretical Physics Computing Facility at Brown 
University and the National Energy Research Sci- 
entific Computing Center. This work is sup- 
ported by DOE contract DE-FG0291ER40688, 
Tasks A and D. 
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